Non-thermal radiation from a runaway massive star 
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Abstract: We present a study of the radio emission from a massive runaway star. The star forms a bow shock 
that is clearly observed in the infrared. We have performed VLA observations under the assumption that the 
reverse shock in the stellar wind might accelerate charged particles up to relativistic energies. Non-thermal 
radio emission of synchrotron origin has been detected, confirming the hypothesis. We have then modeled 
the system and we predict a spectral energy distribution that extends up to 7-rays. Under some simplifying 
l^j assumptions, we find that the intensity at high energies is too low to be detected by current instruments, but the 
future Cherenkov Telescope Array might detect the source. 
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1 Introduction 

CM 

Runaway OB stars (Gies and Bolton 1986) can produce so-called stellar bow shocks on the surround- 
ing interstellar medium. Bow shocks develop as arc-shaped structures, with bows pointing in the same 
direction of the stellar velocity, while the star moves supersonically through the interstellar medium 
(ISM). The stellar and shock-excited radiation heat the dust and gas swept by the bow shock. The 
dust, in turn, re-radiates the energy as mid to far IR flux. 

Van Buren & McCray (1988) looked for bow-shaped features near high- velocity O stars and found 
an IR candidate close to the O supergiant BD+43°3654 (a, <5[J2000] = 20 h 33 m 36.077 s , +43°59'07.40"; 
l,b = 82.41°, +2.33°). Comeron & Pasquali (2007) related the star BD+43°3654 to a bow shock de- 
tected with the Midcourse Space experiment (MSX) at D and E bands. Also, data from the NRAO- 
VLAQNVSS Survey (Condon et al. 1998) revealed a coma-shaped source of ^ 7 arcmin spatially 
coincident with the MSX feature (Figure [TJ. 

A radio study of the bow shock can shed light on the physical processes that give rise to high- 
energy emission from a stellar source, regardless of the history of the runaway star. The shock can 
accelerate particles up to relativistic energies by Fermi mechanism. Energetic electrons will cool 
through synchrotron radiation, producing a non-thermal radio source. We carried out radio observa- 
tions at two frequencies to study the nature of the emission from the bow shock of BD+43°3654. 



National Radio Astronomy Observatory - Very Large Array (http://www.vla.nrao.edu/). 



"The multi-wavelength view of hot, massive stars"; 39 Liege Int. Astroph. Coll., 12-16 July 2010 
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Figure 1: MSX-D band image (color scale) superposed to 1.4 GHz-NVSS contours. Levels are: —2, 
2 (2a), 5, 8, 11, 15, 19, 24, 29, 50, 70, and 90 mJy beam" 1 . 

2 Observations and results 

Our continuum observations were carried out with the Very Lare Array (VLA) at 1 .42 GHz (C config.) 
and at 4.86 GHz (D config.). Figure 2 presents the resulting images after primary beam correction, 
re-gridded with the same synthesized beam of 12". There is emission at both frequencies along the 
extension of the MSX source. The hypothesis of a physical association between the star and the 
radio/IR features is supported by the very good agreement of the residual proper motion of the star 
and the direction from the star to the apsis of the bow shock (Fig. [2]). We used the continuum images 
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Figure 2: Continuum emission at 1.42 GHz (left), and at 4.86 GHz (right). Contour levels are —3, 3, 
6, 10, 15, 20, 25, and 60 times the rms of 0.3 and 0.2 mJy beam" 1 . BD+43° 3654 is marked with a 
cross. The arrow represents the velocity vector of the star, derived from proper motions corrected for 
local motion of the surrounding ISM (see text). Synthesized beams of 12" x 12" are shown in the top 
right corners. 

at 1.42 and 4.86 GHz to build a spectral index distribution map. We only considered input pixels with 
a signal-to-noise ratio > 4. Besides, the spectral index map was masked for a signal-to-noise ratio > 
10. Figures 3 and 4 show the spectral index distribution and corresponding noise maps. 
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Figure 3: Spectral index distribution. 
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Figure 4: Spectral index error distribution. 

3 Bow-shock emission 

Most of the area shows a source of non-thermal radiation with index < a >= —0.5 (S u oc z/*), as 
obtained from the VLA data. We adopted a distance to the bow-shock of 1.4 kpc (Hanson 2003). 
The distance from the star to the bow- shock is R = 5', or 2 pc. The volume occupied by the bow- 
shock is ~ 4.6 pc 3 . We took a particle density of the ISM in the bow shock region of 100 cm -3 (see 
Kobulnicky et al. 2010). 

The non-thermal radiation is expected from synchrotron emission generated by relativistic elec- 
trons accelerated either at the forward shock in the ISM or in the reverse shock in the stellar wind. 
We estimated the particle energy distribution (n) using the observed flux density and spectral slope, 
and assumed equipartition between magnetic and relativistic particles energy density. 

We considered that the energy density of relativistic particles has three contributions: 

u = u ei + u p + u e2 = J E ei n ei (E ei )dE ei + J E p n p {E p )dE p + J E e2 n e2 {E e2 )dE e2) (1) 

where e l5 p, and e 2 stand for relativistic primary electrons, protons, and secondary electron-positron 
pairs (i.e. pairs coming from charged pion decays), respectively. The relation between primary elec- 
trons and protons energy density is u p = au ei , with a > 0. Three cases were considered: a = (just 
electrons), a = 1 (equal energy density in both species), and a = 100 (proton-dominated case, as 
observed in the galactic cosmic rays). The magnetic field resulted B~5x 10 -5 G. 



The maximum energy of the particles was determined by balancing energy gains and losses. The 
loss mechanisms considered were (i) synchrotron radiation, (ii) relativistic Bremsstrahlung, (Hi) par- 
ticle escape from the radiation region due to convection by the stellar wind, and (iv) inverse Compton 
(IC) scattering of IR, stellar and cosmic microwave background photons. In the case of protons, 
the only relevant losses are proton-proton (pp) inelastic collisions and convective escape. Diffusion is 
negligible in comparison to convection in this situation (the respective timescales are i conv ~ 6 x 10 6 s 
and i di ff ~ 10 13 /(£'/erg) s). Both primary electrons and protons reach energies up to ~ 10 13 eV, 
which is imposed by non-radiative losses, except for a = 100, where synchrotron losses dominate for 
electrons. Figure 5 shows the losses for electrons and protons in the case a = 1. 




Figure 5: Left: acceleration ('Acc'), escape ('Esc'), and cooling times for electrons, due to syn- 
chrotron radiation ('Synchr'), IC scattering of dust photons (TC (dust)'), stellar photons (TC (star)'), 
and backgroud photons (TC (bg)'). Cooling time for relativistic Bremsstrahlung radiation is indicated 
as 'Rel B'. The figure is for the case with equal energy density in electrons and protons (a = 1, see 
text). Right: acceleration, escape, and cooling time for protons due to the pp radiation process ('p-p'). 



4 Discussion and perspectives 

The presence of highly relativistic particles in a dense medium with high photon density can result 
in the efficient generation of 7-rays. The corresponding 7-ray emissivity can be calculated using the 
delta-functional approximation (e.g. Aharonian & Atoyan 2000, Kelner et al. 2006). 

In Fig. 6 we show the spectral energy distribution obtained for the case a = 1, with all con- 
tributions included (synchrotron self-Compton losses are negligible). The total luminosity from pp 
interactions is similar to that obtained from relativistic Bremsstrahlung of electrons. The IC up- 
scattering of IR photons is the major contribution at high energies, with a peak around 100 GeV. The 
detectability of the source by the Fermi LAlQ 7-ray observatory will depend on the actual particle 
density and the contribution related to the secondary electrons at large a. The pp contribution extends 
well into the TeV regime, but its weaker and will be difficult to detect with the current ground-based 
Cherenkov telescope arrays. 

For the case a = 100, the relativistic particle contain is proton-dominated and 7-rays from pp 
process dominate the high energy spectrum. The CTA^] North observatory might detect the source 
yielding information on the high enery cutoff. 

2 Fermi Large Area Telescope (http://www-glast.stanford.edu/). 
3 Cherenkov Telescope Array (http://www.cta-observatory.org/). 



Observations of the spectral slope at high energies can be used to identified the proton content 
through the luminosity level, and the proton spectral index. Radio polarization data will provide 
additional information of B. X-ray observations will allow to determine the cutoff of the synchrotron 
spectrum, directly related to the maximum energy of the electron^} This, in turn, would yield valuable 
information on the actual value of B and the correctness of the equipartition hypothesis. 
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Figure 6: Spectral energy distribution for the case a = 1. Acronyms as in Figure 5. Measured radio 
fluxes from VLA observations ('VLA Obs') and MSX luminosity at D band are also represented. The 
contribution from secondary pairs is negligible in this case, so is not shown here. 
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4 Notice that the situation is quite different from that of colliding winds, where the particle acceleration occurs in a 
region of high photon density, with dominance of IC losses. 



